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1. Introduction 

String theory has been successful in resolving some longstanding problems, such as 
the existence of a consistent theory of quantum gravity. However, many problems remain 
unsolved [Q . One of the most important problems is that by far string theory has not made 
any concrete predictions verifiable by experiments, thus we do not know whether string 
theory is a realistic physical theory or not. Cosmology may be an important arena to 
test string theory. For example, any evidence of the existence of topological defects such 
as cosmic superstrings, can be an important support for string theory. Cosmic strings 
can have two different origins, the field theoretic one and fundamental string theory. The 
discovery of fundamental cosmic strings would be a spectacular way to verify string theory 
||. Cosmic strings from string theory are characterized by some properties not shared by 
GUT cosmic strings [[|. 

In 1980's, it was generally believed that the perturbative fundamental strings can 
not become cosmic strings, due to an argument of Witten [|J. Moreover, the tension of a 
fundamental string is close to the Planck scale, while cosmic strings with such a tension are 
ruled out by experiments. Recently, research of compactifications in string theory shows 
that the string tension measured in the four-dimensional Einstein frame can be much 
smaller, and in some situations the instability problem is evaded. Thus cosmic strings as 
fundamental strings may indeed exist and can be observed in the future experiments. 

The current research interest of the creation of strings at the end of inflation is focused 
on the investigation of the final results of collision of branes. Cosmic strings are inevitably 
produced in this process as topological defects. Research along this direction is spearheaded 
by Polchinski and Tye, and their collaborators || @ ■ 

In this paper, we will study creation of strings in a more traditional fashion, namely, 
through gravitational pair production in a time-dependent background. Related work has 
been done by Gubser with an effective field theory viewpoint, and a steepest descent 
contour method has been developed to estimate the production rate of the strings. 

As discussed in 0, in a regime of parameter space where a spacetime description 
gives a good approximation of string dynamics, the on-shell constraint for a given string 
state boils down to a differential equation describing an oscillator with a time-dependent 
frequency. When the quantization of strings is carried out in a particular background, the 
"frequency" uj{t) is determined. The quantization of strings in a de Sitter background was 
recently done by Li et al. ||, and the spectrum of such "small strings" is obtained. We 
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will use the method developed in to estimate the total rate of string creation, using 
the spectrum i.e. the equation of motion of string state obtained in ||, which is different 
from that of . The original equation of motion derived in || is in de Sitter space where 
the Hubble parameter is a constant, but as we will show that in fact it is also valid at the 
end of inflation and during reheating, where the Hubble parameter is a function of time 
instead of a constant (see Appendix A). We will show this in Appendix B. 

The main result of our investigation can be summarized as follows. Strings are gener- 
ally produced gravitationally at the end of inflation and during reheating, and the energy 
density of strings produced is highly suppressed by an exponential factor multiplied by an 
power factor. Since the Hubble parameter is much smaller than 1 in the unit of a' = 1, 
this energy density is very small. Our estimate is approximate quantitatively, due to the 
fact that there are some approximations used in deriving the string spectrum derived in 
U and in our analytic method in estimating the string production rate. However, this 
semi-quantitative result strongly suggests to us the picture that highly excited strings are 
hardly produced during reheating and the production rate is very small. 

The organization of this paper is as follows. In section 2 the density of string states 
is calculated, based on the spectrum obtained in [Q. In section 3 the creation rate and 
the energy density of strings are estimated. The final section is devoted to discussions. In 
Appendix A we solve the Friedmann equation directly to get the Hubble parameter as a 
function of time at the end of inflation and during reheating. In Appendix B We show 
that the spectrum and equation of motion of strings obtained in f|] are valid in a general 
flat FRW background, not only in pure de Sitter space. 



2. Degeneracy of String States 

The spectrum of strings in a near de Sitter background is different from the one in 
flat spacetime, and depends on two integers (we shall consider the bosonic sector only in 
this paper), as will be shown in the next section. These integers are eigenvalues of two 
operators: the number operator N and the other operator L, defined respectively as follows 

d oo 

i=lm=l ' . . 

d oo V ' 

i=\ m=l 
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where d = D — 1 and D is the number of dimensions of spacetime, and are the 
occupation numbers of the left- mover and the right- mover respectively, i is the space index 
and m is the oscillator index. We denote n as the eigenvalue of N and / as the eigenvalue 
of L. The degeneracy of states at level n and fixed I is denoted by -D n ,z, which is encoded 
in a generating function as the coefficient of z n w l 

Z(z,w) =trz N w L , (2.2) 



or 



Z(z,w)= D n,iz n w l . (2.3) 



,1=0 



Using ( |2.1| ) and (|2.2| ) , the generating function can be evaluated by an elementary method 



of quantum statistical mechanics as follows 

d oo oo 



Z(z, w) = J] fj E z m ^ +N ^w N -+^ +2N ^. (2.4) 

i=l m=l Kfi fji _n 



Summiner over 



Z(z, w) 
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z mNl 



(2.5) 



nn e rr^/c^)' 



i=l m=l TV* =0 



where 



(1 - z m )w N ™ 

f(w,Nl) = ± - J . (2.6) 
1 — ^ m w /iV «i +i 



If w = 1, (|2.5|) is just the usual generating function of the degeneracy of bosonic string 
states [H . The effect of w is to deform this generating function. This deformation is small 
when it comes to evaluate the coefficient D n ^ by the steepest descent method. We will 
expand the generating function i.e. ( |2.6| ) near w = 1 to the second order, with higher 
order terms truncated. This assumption is reasonable as we will show in the saddle point 
calculation. 



f(w,Ni) = l + (w-l) 



l J rZ m 

+ N l ' ' 

I - Ly 7J 



l-z m " L \l-z 



+ 0((w-lf). (2.7) 



Of course one can sum over NL first. 
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Thus the summation in (|2.5| ) reads: 

E T^ f{w ' N ™ 



m 

then the generating function is 



(1 - z m ) 2 



\ + {w-\) 



2z r 



(1 - z m ) 2 



(2.1 



n (T 

m=l v 



1 



,m\2 



l + (w-l) 



2z J 



(1 - z m ) 2 



(2.9) 



Let z = exp(— ^) and x = Sr. Thus x is continuous when T is large enough. In 
our following calculation, around the saddle points, z is close to 1 thus T is indeed large. 
The infinite product in ( |2.9| ) is approximated by an exponential of an integral. Taking 
logarithm of both sides of (|2.9|) we obtain 



\nZ(z, w) 
d 



-2T 

7T 2 T 
7T^ 



dx ln(l - e~ x ) +T dxln 1 + (w - 1 



2e" 



(1 



-aA2 



+ 2(w - 1)T 



dxe~ 



i (1 



(2.10) 



2(w-l)z 



3 In 2 (1 — 2;) In z 

in the second line we have truncated higher orders of w — 1 in the logarithm function. 
By definition, the degeneracy of string states is 



D T 



dw 
2tH 



dz Z(z, w) 



27tiz n+1 w l+1 . < " 2 ' 1 ^ 
Here Z(z,w) vanishes rapidly as z — > 1 when w < 1, and z n+1 is very small for z < 1 
when n is very large. Consequently, for large n, there is a sharply defined saddle point for 
z near 1. Indeed, the factor i 



is stationary for 



exp[ 



lnz 



2tt 2 A(w-1)z 
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(n + l)lnz - (I + l)\nw] 



1 



n+ 1 

Therefore one finds that as n — > 00 



2tp 



lnz 



1 2?r 
n+ 1~3~ 



D n 1 ~ 71 5/4 exp 



2tt 



6nJ 5 ( Z 



6n 



7T 



(2.12) 



(2.13) 



(2.14) 



where the 5- function comes from the integration over w, and the exponential factor is the 
ordinary degeneracy of bosonic string states in four-dimensional spacetime. 



Here d = 2 is the dimension of physical states, i.e. the transverse oscillators are iV* 



1,2. 
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3. Energy Density of Cosmic Strings 

In this section we will use the steepest descent contour method developed in J7| [|10| to 
estimate the energy density of strings produced during reheating. In ||, the quantization 
of bosonic strings has been done in a de Sitter background. The on-shell constraint for 
quantum states of a string leads to an equation of the form: 

jd t 2 + 3Hd t + k 2 e~ 2Ht + 4N + 2E 
- J2 H 2 (l + 2N* m N* m + + jty }<p(N^ A> , u,, V) = 0, (3 ' 1} 

m,i 

where we take a' = 1. H is the Hubble parameter, N^, are occupation number 
operators and N = m(N^ n + NJ n ), k % is the momentum vector in the four-dimensional 
spacetime and Eq is the center-of-mass energy. A general physical state \<j>) corresponding 
to the string modes is related to </>(A^, N^, u>, k l ) as follows, 

10) = E (3.2) 

AT* ,N* 

where the definition of | N^, N^, u, k l ) can be found in ||. Eq.fl3.lp is regarded as the 
equation of motion for the field of the corresponding string state. As we mentioned in the 
introduction, this equation is different from that used in , this difference makes our new 
result different from others. 

There is no string production in pure de Sitter space with a constant Hubble param- 
eter, even with the modified string spectrum as in fl3.1|) . The only chance for string pro- 
duction to occur is the short period of reheating during which H becomes time-dependent. 
Thus we have to make a step forward, i.e. to generalize the equation of motion (|3.1| ) to the 
case when H varies with time. This is developed in Appendix B. From now on {k 1 ) 2 e~ 2Ht 
in ( p.lp is replaced by (k l ) 2 /a(t) 2 , where a(t) is the cosmological scale factor. 

It is convenient to introduce 4>{t) via $(£) = a(f)(t), thus the equation of motion for 
<f>(t) is: 

where 

C2 = E ( 3 " 6 + N ™ + + 2N ^m) (3-4) 

m,i 



a(t) 



+ AN + 2E - C 2 H< 



= 



(3.3) 
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and 

e = - w (3.5) 

denned as the so-called slow-roll parameter, which is roughly equal to 1 at the end of 
inflation, and the dot denotes the derivative with respect to t. The slow-roll parameter is 
not larger than 2 in the case under study, see Appendix A. 

In conformal time rj defined by adr\ = dt 1 we can eliminate the first order derivative 
term, and we will use prime to denote the derivative with respect to the conformal time. 
Thus (|3.3| ) becomes: 



4>{ji)" + W(ri) 2 (j>(v) = 0, (3.6) 

where 

W{r]) 2 = k 2 + (4N + 2E )a(r]) 2 - C 2 a(r]) 2 H (r?) 2 . (3.7) 

Having obtained the equation of motion ( |3.6| ) , now we use the steepest descent method 
to extract the approximate string pair production rate from (|3.6| ) . The steepest descent 
method was developed by various authors, especially, it was used to estimate the string 
production rate by Gubser [[7] . The key assumption is that the occupation number \(3\ 2 
for a given mode is always much less than 1, where (3 is the Bogliubov coefficient. Setting 

m = afa) -j /" du W(u) m j P du W(u) m 

with the requirement |ct(r7)| 2 — \P(rj)\ 2 = 1, we recast the equation ( |3.6|) into 

// s W 2i P duWiu) Q ( \ alt \ W' -2i f v duW(u) / \ , Q n s 

Using the assumption [3(i]) 1 and 01(77) ~ 1, we obtain an approximate formula for (3 

/3^y dr] — exp{-2i J duW(u)J . (3.10) 
The integral in the exponential of ( |3.10| ) can be calculated as follows 



duW(u) 



Vi 



v rv 

+ . 

Vi J v 



duW(u) (3.11) 



where rji is some initial time and 77* is denned to make W(rj*) = 0. Here the second term 
on the right hand side of (|3.11|) can be calculated as follows: 



/ duW{u) = / du^k 2 + (AN + 2E )a(r]) 2 - C 2 a{r)) 2 H(r)) 2 

J rj* J r)* 

w / du^J[{AN + 2E Q )2aa' - C 2 2aa'H 2 - C 2 a 2 2HH'](u - rj*~) (3.12) 

= -5 3/ V( 4Ar + 2E )2aa' - C 2 2aa'H 2 - C 2 a 2 2HH' 
3 

where we have expanded terms in the square root around rj*, and defined 6 = 77 — 77*. Thus 

/3^/ exp(-2z f duW(u)) (3.13) 



where 



and 



// ?? ^ eXp( "^ 3/2 ^ ) (3 ' 14) 



S(r)*) = (4N + 2E Q )2aa' - C 2 2aa'H 2 - C 2 a 2 2HH'. (3.15) 
Expanding W(rj) around 77*, we get 

/o= 4 /lf exp H 53/2vWy )- (3 - i6) 



From |T0[ we know that integrals such as ( |3.16[ ) can be calculated by the contour integral 
method, and the result for (|3.16| ) is simply I = in/3. 
Now ( |3.13|) can be written as: 



3 



r \ I rrj 



(3 w — exp ( -2i I duW{u) ) exp I -2% I duW{u) ) (3.17) 



Vi / V J r 



where r is the real part of 77* = r — ut, and r, u are real with u > . Since what we need 
is the modulus of f3, the main contribution comes from the second exponential function in 
( ^.17| ) whose argument is the following integral and can be expanded as: 

r v * (-in) 2 (-iv) 3 

/ duW{u) = W(r)(-iu) + W\r)^-^- + W"(r)^— !- + ■■■. (3.18) 

Jr 2 6 



As long as |^r| <C |^|, we can truncate this expansion to the first term (even terms 
do not contribute to the modulus of (3 since they are real), and indeed in the following 
calculation one will see that | | *C | A | is satisfied in our case. Thus we get 

\(3\ 2 « (|j exp (-4uW(r)) . (3.19) 



Here the imaginary part of ?y* can be derived by using ( |3.7| ) and expanding W(rj*) around 
r as follows! 

=W 2 {rf) = W 2 (r) + 2W(r)W'(r)(-iu) + [W(r)W"(r) + W' 2 (r)] (-iu) 2 

(-iu) 3 ( 3 ' 2 °) 

+ [W'"(r)W(r) + W'(r)W"(r) + 2W'(r)W"(r)] 

3 

We solve these equations as follows 

W 2 [r) - [W(r)W"(r) + W' 2 (r)]u 2 = 
6W(r)W'(r) - [W" (r)W(r) + 3W" (r)W' (r)]u 2 = o" 

Then 



y/w"(r)/W(r) + (W'(r)/W(r)) 



In our case W(rj) is expressed in (|3.7|) , thus we get 



and 



(3.21) 



j) 2 exp( ; ^ | . (3,2) 



W a V a a a I a 



W \W J \w 
Thus we write approximately: 

L9 fc n, 2 w exp ^ V / r r ; S , 3.25 

where we wrote (3 as function of comoving momentum k and excitation modes n, I explicitly, 
where n and / are eigenvalues of operators in ( |2.1| ) respectively. We have dropped the 



We truncate the expansion up to the fourth term because one can check that the next order 



terms are much smaller than these terms. And one can see that the condition I 1 <C I -4r I is a ls° 
satisfied. 
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factor of (7r/3) 2 , H r and R r correspond to the Hubble expansion rate and Ricci scalar for 
the metric ds 2 = a(n) 2 (drj 2 — dx l ) respectively. We emphasize that all time-dependent 
quantities in (|3.25| ) are evaluated at 77 = r, where r is the real part of 77* given by W(n*) = 0. 
From ( |3.25| ) we can see that indeed the production of highly-excited strings i.e. strings 



with large n and I are exponentially suppressed. 

The total energy density of strings produced may be written as 

= 77^3 I dkk 2 Y, D n,i \Pk(nJ)\ 2 M^fa), (3.26) 
2n 2 a{r,) J nJ 

where k is the comoving momentum, M n j (n) is the energy of a single string with excitation 
modes (n,Z), given by 

M n M 2 = 4iV + 2E -J2 H2 (l + + < + N^) 

m,i 

= An + 2E -H 2 (l + l) ■ ( 3 - 27 ) 

ss An - H 2 l 



Now insert ( |3.25|) and (|3.27| ) into ( |3.26|) , we obtain approximate formulas: 



P ^ ^ 27To3 / dkk2 ^l2 n * eXP (y V/ ^) ^ ~ 7T^ X 

/ k 2 / a 2 r + An - a 2 r H 2 l\ r — 

exp -4 — ' r r r \^An - H 2 l 

1 \ y/4n(H* + R r /6) J 

1 [„,2V^ -5 /27T , \ 

~ — / dkk n 4 ex P I V6n J x 

n \ / 



(3.28) 



^An{H 2 + R r /6) J V 7T 2 



For highly excited string states, neglecting A; is a good approximation, though not a uniform 
one if a{rj) becomes arbitrary small in the past. We get: 

f\ 1 f,n2 C max a -a / f8-12a 2 r H 2 /7r 2 2V6n\ \ 

p(v) ~ / d/c/r / dnn 4 exp <^ - — r n = vn > . (3.29) 

Py,, a*J A P \ \ ^H 2 + R r /6 3 / J 

where we have dropped the constant factor. We define A = ( 8 j^g r H r Jn^ — 2yW ] f or 
short, the above integral can be approximated as: 

p(??) ~ ^4 / ^ 2 ( Er / (v^iL*) - £r/ (v 7 !)) , (3.30) 
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where Erf(x) is the error function. For very large x, approximately we have 



Erf{x) ss 1 - 



X\ TV 



1 + f ( _ ir l-3-5-(2n-l 



77 = 1 



X\ 7T 



(3.31) 



Thus we have 



where A is given by A 



a(r)) 3 A 3 /< 



-Ay/r 



S-12a 2 r H 2 r /-K 2 



e — 



1/4 
Tlmax 



(3.32) 



, , , . Since we consider the highly excited 

y/H?+R r /6 3 J 6 J 

strings, i.e. strings with small momentum k, the integral in ( |3.30| ) with respect to the 
comoving momentum k will contribute a small factor which we have dropped. This result 
will not affect the qualitative behavior of the production rate of strings with respect to the 
string excitation modes n and I. The upper limit of n is roughly n max ~ H~ 2 /4 ^> l(See 
Appendix B), thus p ~ , p^ 3/2 . We emphasize that indeed H r <C ^==. Thus in unit 
where a' = 1, if r 1 and A 3> 1, the above approximation is qualitatively correct. 
Especially, since A>1, from ( |3.32| ) we can see that the energy density of strings produced 
is very small and exponentially suppressed, i.e., highly excited strings are hardly produced 



m our case. 



4. Discussion 

We have estimated the energy density of strings produced at the end of inflation and 
during reheating, our main result is 

here A is given by A = ( 8 12a r H r a 1^ _ 2^2L ] \Y e have reinstated a' which has been 

yy/H^a'+Rra'/6 J J 

set to 1 in our paper. Although it is difficult to get the explicit form of H r due to the 
complicated equations of ( |3.21|) , one can make sure that H r \foi' must be much smaller 
than 1, i.e. the Hubble scale in string production process is much lower than the string 
scale. In other words, the curvature radius is much larger than the string length, and a 
spacetime description gives a good approximation. In the case of a small H r , due to the 
large exponential factor in ( |f .1|) , one can see that the energy density of strings produced 
is exponentially suppressed and indeed highly excited strings are hardly produced. (|3.7|). 
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Planck scale, otherwise the effective field theory viewpoint we used will broken down, so 
that the energy density is also small. 

In conclusion, we have shown in this paper that highly excited strings are hardly 
produced at the end of inflation, because \(3\ 2 is highly suppressed by a exponential factor 
within and the degeneracy of highly excited strings is not sufficiently large to compensate 
it, thus the energy density is also suppressed by this exponential factor. 
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Appendix A. Hubble Parameter During Reheating 

In most of popular inflation scenarios, the temperature is practically zero during 
inflation, relativistic matter is produced during the short reheating period when the inflaton 
oscillates coherently and decays to matter. Generally it is not known how the inflaton is 
coupled to a generic string state, so the usual reheating mechanism is not easily applied 
to the production of strings. 

However, the spacetime metric is also coupled to strings, the details of the coupling 
can be seen from the string spectrum directly. When H remains nearly a constant, there 
is no string production. During the reheating period, the Hubble parameter is no longer 
a constant, and can be estimated by solving the Friedmann equation. In solving this 
equation, we should also take radiation into account. A more rigorous treatment should 
also take strings produced in the process into account, however, we do not know how to 
compute string energy density as a function of time (to this end, it is required to compute 
the string production rate per unit time). The Friedmann equation reads 



where we set 9,-kG = 1, and p<j>(t) and Pj(t) are the energy densities of the inflaton and 
radiation respectively, their equations of motion are 



3H 2 = p (f> + p 7 , 



(A.l) 



+ ZHp^ + Tp^ = 0, 



(A.2) 



and 



p 7 + AHp 7 - Tptj, = 0, 



(A.3) 
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where V is the decay rate of the inflaton and dot denotes the derivative with respect to 
the comoving time t. 

Taking derivative of ( |A.1| ) with respect to t and using ( |A.2| ) ( |A.3| ) to eliminate p^ and 



p 7 , and then using ( |A.1| ) again, we find 



2e " 3 ,K AS 

Pi = J^eP^ ( A - 4 ) 
where e = —H/H 2 . Now the Friedmann equation is 

3H 2 = j^- e P4>- (A.5) 

Combining ( |A.4j ) and ( |A.5| ), we get 1.5 < e < 2. This result can be generalized when there 
are more energy components in the Universe. 

We solve (|A.2| ) in the limit@ of T ^> H, which means that the inflaton decays very 
fast, indeed this is the case during reheating. Thus (|A.2|) is simplified to 



P<t> + T P(j} = 0. (A.6) 

The general solution is 

p^ = p exp (-Ft) , (A.7) 

where po is an integration constant. Using solution ( |A.7| ), ( |A.5| ) can be rewritten as a 
differential equation of H explicitly 

6H + 12H 2 - p exp(-Tt) = 0. (A.8) 

This equation can be cast into the standard form of Bessel equation by changing variable 
t into r = exp(— Tt/2), 

The general solution to this equation is a linear combination of the modified Bessel func- 
tions of the first kind Iq and of the second kind Kq: 

" 2 = c ^(^t) +c ^{^t)' (A10) 

where c± and C2 are integration constants. 



Of course one can exactly solve the equation, but that is not necessary. 
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Appendix B. Constraint Equation with Arbitrary H(t) 

The equation of motion of bosonic string states derived in || is valid in pure de Sitter 
space, of which the Hubble parameter is a constant. In order to study string production, 
we focus on the reheating phase, where the Hubble parameter is varying with time, the 
spectrum formula of || cannot be used directly for our purpose. In this appendix we 
generalize the original result of |J to the case of arbitrary H(t) = a(t)/a(t). We refer the 
readers to the original paper H for more details of deriving the equation of motion when 
H is a constant. 

The key point of the generalization is to re-calculate c m in (3.16) of || in arbitrary 
H(t) case. The general definition of c m is 



-,i{4>m—'>Pm,) 



*2 -2 

a e 



/t ft 
duXm^) _ p*2^+2i J du\ 7n (u) 



x r 



2Ar 



(B.l) 



where the most general form of a m and (3 m are given by (3.7) and (3.8) of 

a m = cosh( 7m )e l,5 ™+ 1 ^ , a m = cosh( 7m )e l<5 ™+ 1 ^ , 



P m = sinh(7 m )e l 
X m in QB.lj ) is defined as 



(3 m = sinh(7 m )e 



(B.2) 



A m = sgn(m)\ - r]dfr] 



(B.3) 



where rj = eH t^ as in ||, and should be replaced by a formula in which exp(Ht) is 
replaced by a(t), and 7 m and 5 m in ( |B.2| ) can be solved directly from (3.9-3.10) of ||. 



Jr. 



cosh 1 



AmXr 



2 2 

1 + A m H 2 



1 

+ 2 



(B.4) 



arctan 



— 2A TO .r 



du\ m (u). 



(B.5) 



2r 2 -r. 

here we define r = 77/77 for short. Insert ( |B.2| ) into ( |B.1| ) and use ( |B.3| )-( PT5|) , after a 
tedious calculation we finally get 



m ) I 



t + M£)-(t) 2 l-^(! + £) 



a 1 a I di \ I a\ i 4m 2 fa, 6j \ 
a + 0t \2^) ~ [a) + —{a + 2^) 



(B.6) 
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What is needed in the calculation of the main context is the modulus of c r 



\Cm\ = ~ = H(t), 



(B.7) 



it agrees with the original result of || where H is a constant. (|B.7|) tells us that the equation 
of motion of string states in the reheating phase where Hubble parameter is varying with 
time H = H(t) is simply given by the equation of motion of |§ when H is replaced by 
H{t) = a(t)/a(t). 

At the end of this appendix we want to recall that the real condition of each A m is 



__ (ff+ _ )2 _ at( _ )>0 



(B.8) 



where uj 2 = 4N + 2Eq + (p l /a) 2 . In fact if Ai is real, so are X m for m > 1. Then the 
condition becomes 



1 



Because 



(V_ 



2 \au J 4 \auj 



AN + 2E + (p i /a) 2 
and e G (1.5, 2), the condition ( ]B.9| ) is approximated by 



H 1 > 0. 



< 1 



(B.9) 



;b.io) 



lo- 



l+ e ~ 

2 \auj 



H z > 0. 



Thus we get the upper limit of n as 



n < n max = - (h~ 2 - (l + 1) k 2 - 2E 



:b.iv, 



(B.12) 



where k l = p l /a is the physical momentum. 



14 



References 



[1] J. Schwarz, |hep-th/07022T9 

[2] 

[3] 



[4] 
[5] 
[6] 
[7] 
[8] 
[9] 

[10] 



J. Polchinski, AIP Conf. Proc. 743, 331 (2005) [Int. J. Mod. Phys. A 20, 3413 (2005)], 
|hep-th/0410082 . 



A. Vilenkin and E. Shellard, Cosmic Strings and Other Topological Defects, Cambridge 
University Press (Cambridge, 2000); M.Hindmarsh and T. Kibble, Rep. Prog. Phys. 
58, 477 (1995). 

E. Witten, Phys. Lett. B153, 243 (1985). 

N. Jones, H. Stoica, and S. Tye, JHEP 07, 051 (2002), |hep-th/0203163 . 

S. Sarangi and S. Tye, Phys. Lett. B536, 185 (2002), |Kep!h/0204074j . 

S. Gubser, Phys.Rev. D69 (2004) 123507, |hep-th/0305099| ; |hep-th/0312321 . 



M. Li, W. Song and Y. Song, |hep-th/0701258| , JHEP 04 (2007) 042. 

M. Green, J. Schwarz and E. Witten, Superstring Theory, Vols. I and II, Cambridge 

Univ. Press (1987). 

E. Brezin and C. Itzykson, Phys. Rev. D2 (1970) 1191; D. Chung, |hep-ph/ 9809489 



15 



